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The  solar  cell  operation  temperature  is  an  important  input  in  models  for  the  electrical  characterization  of 
high  concentration  photovoltaic  cells  and  modules.  However,  the  direct  measurement  of  this  tempera¬ 
ture  is  difficult  in  these  kinds  of  devices.  Because  of  this,  in  recent  years,  the  scientific  community  has 
proposed  different  methods  for  indirectly  calculating  the  cell  temperature  in  high  concentration 
photovoltaic  modules  from  atmospheric  parameters  and/or  easily  measurable  parameters  on  the 
module.  In  this  paper,  a  comprehensive  review  of  existing  methods  for  cell  temperature  calculation  in 
high  concentration  photovoltaic  modules  for  electrical  characterization  is  presented.  The  different 
methods  are  summarized  and  a  comparative  analysis  is  done.  Required  inputs,  advantages  and  technical 
difficulties  of  each  method  are  highlighted.  Also,  an  experimental  campaign  has  been  carried  out  at  Jaen, 
south  of  Spain,  in  order  to  quantify  the  accuracy  of  the  methods.  Results  show  that  methods  based  on 
direct  measurements  on  the  module  are  somewhat  more  accurate  than  methods  only  based  on 
atmospheric  parameters.  The  choice  of  the  most  suitable  method  for  a  specific  application  will  depend 
on  the  availability  of  module  information,  on  the  required  accuracy  and  on  technical  issues. 
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1.  Introduction 

High  concentration  photovoltaic  (HCPV)  systems  use  optical 
devices  (lenses  or  mirrors)  to  concentrate  the  sunlight  onto 
small  solar  cells.  Although  there  are  a  lot  of  possible  configura¬ 
tions  in  order  to  implement  this  concept  [1,2],  a  typical  HCPV 
system  consists  of  a  two-axis  solar  tracker  on  which  HCPV 
modules  are  mounted,  each  module  composed  of  several  elec¬ 
trically  connected  solar  cells  with  their  associated  optics.  The 
small  size  of  the  solar  cells  allows  incorporating  high  efficiency 
III— V  multi-junction  solar  cells,  while  other  silicon-based  solar 
cells  can  also  be  used  [3],  One  of  the  advantages  of  these  systems 
is  the  high  potential  of  costs  reduction  because  they  replace 
expensive  semiconductor  materials  by  cheaper  optical  materials 
[4],  As  well,  the  operation  with  high  efficiencies  implies  an 
increase  of  produced  energy  for  the  same  installed  area  than  a 
conventional  photovoltaic  system  [5], 

The  operation  temperature  of  the  cells  in  an  HCPV  module  affects 
its  performance.  This  is  because  the  band-gap  of  the  semiconductors 
that  compose  the  solar  cells  reduces  when  temperature  increases.  This 
implies  a  change  in  the  solar  cells  electrical  parameters:  open-circuit 
voltage,  fill  factor,  maximum  power  and  efficiency  decrease  with 
increasing  temperature,  while  short-circuit  current  increases  (Fig.  1). 
These  effects  have  been  observed  for  both  silicon-based  solar  cells  and 
III— V  multi-junction  solar  cells  [6-9],  As  HCPV  modules  behavior  is 
influenced  by  the  cells  inside,  HCPV  modules  electrical  parameters  are 
affected  by  temperature  in  the  same  way  [10,11],  Thus,  knowledge  of 
the  cell  temperature  is  critical  for  characterizing  the  behavior  of  HCPV 
modules. 

While  cell  temperature  is  an  important  input  in  models  for  the 
electrical  characterization  of  HCPV  cells  and  modules  [12],  the  direct 
measurement  of  this  temperature  in  HCPV  modules  seems  to  be 
difficult  because  it  requires  accessing  inside  the  module  and  placing  of 
a  small  temperature  sensor  very  close  to  the  solar  cell  [13].  Although  in 
conventional  flat-plate  modules  the  cell  temperature  can  be  ade¬ 
quately  estimated  by  measuring  the  temperature  on  the  back  of  the 
module,  in  HCPV  modules  the  complex  thermal  behavior  and  the  use 
of  heat  exchangers  for  keeping  the  cell  temperature  in  acceptable 
levels  and  avoiding  degradation  make  difficult  this  practice  [14,15], 
Because  of  these  difficulties,  in  recent  years  the  scientific  community 
has  devoted  efforts  in  developing  methods  for  indirectly  calculating 
the  cell  temperature  in  HCPV  modules  from  different  parameters 
(atmospheric  parameters  and/or  easily  measurable  parameters  on  the 
module). 

The  cell  temperature  is  not  only  important  in  the  modeling  of 
HCPV  modules,  but  it  has  interest  in  the  power  rating  of  these 
devices.  Conventional  photovoltaic  modules  are  rated  at  specific 
conditions  of  cell  temperature,  usually  at  25  C  of  cell  temperature 
[16].  In  contrast,  HCPV  modules  are  usually  rated  at  defined  condi¬ 
tions  of  environmental  parameters  due  to  the  difficulty  of  measuring 
the  cell  temperature.  For  instance,  the  ASTM  standard  E2527-06  [17] 
requires  ambient  temperature  of  20  °C,  direct  normal  irradiance  of 
850  W/m2  and  wind  speed  of  4  m/s.  Although  establishing  these 
conditions  avoids  the  need  of  measuring  cell  temperature,  the 
obtained  rated  maximum  power  is  not  directly  comparable  to  that 
obtained  for  other  photovoltaic  technologies.  This  is  an  issue  for 
market  penetration  of  HCPV  technology  [18],  Thus,  it  would  be 


desirable  to  have  rating  procedures  based  on  standard  cell  tempera¬ 
ture  conditions  [19]  and  methods  that  reliably  calculate  cell 
temperature. 

In  this  paper,  a  comprehensive  review  of  existing  methods  for 
calculating  cell  temperature  in  HCPV  modules  for  electrical  char¬ 
acterization  is  presented.  Its  aim  is  to  help  the  photovoltaic 
researchers  and  professionals  in  the  choice  of  the  most  suitable 
method  for  each  application.  The  review  comprises  all  the  relevant 
contributions  in  this  field,  mainly  developed  in  the  last  3-4  years. 
Two  remarks  must  be  done: 

There  are  a  lot  of  methods  for  calculating  the  cell  temperature 
in  conventional  flat-plate  modules  [20],  Among  these  methods, 
only  methods  that  have  been  applied  for  HCPV  modules  are 
included  in  the  present  review. 

In  recent  years,  network  thermal  models  have  been  proposed 
that  predict  temperature  distributions  across  the  module  [21-25], 
These  models  are  oriented  to  thermal  management  and  module 
design  studies.  As  their  implementation  is  complex,  they  are  not 
intended  for  electrical  characterization,  so  that  they  will  not  be 
included  in  the  present  review. 

Taking  this  into  account,  the  reviewed  methods  have  been  grouped 
according  to  the  following  classification:  on  the  one  hand,  methods 
based  on  direct  measurements  on  the  module;  on  the  other  hand, 
methods  based  on  atmospheric  parameters.  Furthermore,  methods 
based  on  direct  measurements  have  been  divided  into  methods  based 
on  module  heat-sink  temperature  and  methods  based  on  electrical 
parameters  (Fig.  2).  Methods  based  on  direct  measurements  are 
expected  to  be  more  accurate  than  methods  based  on  atmospheric 
parameters,  while  methods  based  on  atmospheric  parameters  present 
the  advantage  that  the  cell  temperature  can  be  estimated  at  any 
location  from  meteorological  data. 

The  paper  is  organized  as  follows:  Section  2  is  an  overview  of 
HCPV  technologies;  Sections  3,  4  and  5  summarize  the  methods 
for  cell  temperature  calculation  in  HCPV  modules  according  to  the 
proposed  classification:  methods  based  on  module  heat-sink 
temperature  (Section  3),  methods  based  on  electrical  parameters 
(Section  4)  and  methods  based  on  atmospheric  parameters 
(Section  5);  a  comparative  analysis  of  the  methods  from  both  a 
qualitative  and  a  quantitative  point  of  view  is  done  in  Section  6; 
finally,  Section  7  presents  the  conclusions  of  the  work. 

It  is  important  to  remark  that  not  every  cell  in  an  HCPV  module 
exactly  operates  at  the  same  temperature.  Measurements  on  the  back 
plates  of  HCPV  modules  have  shown  that  cells  near  the  frame  have 
lower  temperatures  depending  on  wind  speed  and  direction  and 
orientation  of  the  module  [26,27],  Because  of  this,  authors  use  to  work 
with  the  “average  cell  temperature"  in  the  module.  This  parameter 
will  be  simply  referred  in  this  paper  as  the  “cell  temperature”  of  the 
module  (Tc). 

2.  HCPV  technology  overview 

HCPV  systems  operate  under  light  concentrations  between  300 
and  2000  suns  [5],  These  systems  can  be  built  with  the  participation  of 
different  conventional  industrial  sectors,  such  as  glass,  steel,  aluminum 
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Nomenclature 

c 

parameter  of  the  method  of  Hornung  (dimensionless) 

Ci,  c2,  c3 

parameters  of  the  method  of  Fernandez  (V/W  m~2, 

Acronyms 

V/°C,  V) 

Cg 

geometric  concentration  (dimensionless) 

ANN 

artificial  neural  network 

d  i,  d2 

parameters  of  the  method  of  Almonacid  (°C/W  m~2, 

ASTM 

American  Society  for  Testing  and  Materials 

°C/m  s_1) 

CEAEMA  Centro  de  Estudios  Avanzados  en  Energia  y  Medio 

DM 

direct  normal  irradiance  (W/m2) 

Ambiente 

DNI* 

reference  direct  normal  irradiance  (W/m2) 

HCPV 

high  concentration  photovoltaic 

he 

short-circuit  current  (A) 

1EC 

International  Electrotechnical  Commission 

I  * 

short-circuit  current  at  reference  conditions  (A) 

1ES-UPM  Instituto  de  ENERGIA  Solar  -  Universidad  Politecnica 

k 

Boltzmann  constant  (m2  kg/s2  K) 

de  Madrid 

L 

thickness  of  a  material  (m) 

ISFOC 

Instituto  de  Sistemas  Fotovoltaicos  de  Concentration 

m 

parameter  of  the  method  of  Hornung  (°C/W  m~2) 

MPP 

maximum  power  point 

n 

effective  diode  ideality  factor  of  an  HCPV  module 

NREL 

National  Renewable  Energy  Laboratory 

(dimensionless) 

PMMA 

poly( methyl  methacrylate) 

diode  ideality  factor  of  a  solar  cell  (dimensionless) 

PVUSA 

photovoltaics  for  utility  scale  applications 

Ns 

number  of  cells  in  series  of  an  HCPV  module 

RMSE 

root  mean  square  error 

(dimensionless) 

P 

electrical  power  harvested  by  an  HCPV  module  (W) 

List  of  symbols 

Pabs 

luminous  power  absorbed  by  the  cells  in  an  HCPV 

module  (W) 

fivoc 

temperature  coefficient  of  the  open-circuit 

P measAVG 

averaged  electrical  power  harvested  by  an  HCPV 

voltage  (V/°C) 

module  while  measuring  an  l-V  curve  (W) 

fivoc.cell 

temperature  coefficient  of  the  open-circuit  voltage  for 

Pmpp 

electrical  power  harvested  by  an  HCPV  module  work- 

a  solar  cell  (V/°C) 

ing  at  maximum  power  point  (W) 

fivoc,cell{X 

c=  1 )  temperature  coefficient  of  the  open-circuit  vol- 

q 

electron  charge  (C) 

tage  for  a  solar  cell  at  one  sun  (V/°C) 

Q 

rate  of  heat  transfer  (W) 

fivocO 

temperature  coefficient  of  the  open-circuit  voltage  at 

Qmeas 

rate  of  heat  transfer  for  a  module  at  open-circuit  (W) 

the  reference  irradiance  (V/°C) 

Qmpp 

rate  of  heat  transfer  for  a  module  working  at  max- 

AT 

temperature  difference  between  the  cell  and  the  heat- 

imum  power  point  (W) 

sink  at  the  reference  irradiance  (°C) 

R 

thermal  resistance  between  the  cell  and  the  heat-sink 

A  Tmeas 

temperature  difference  between  the  cell  at  open- 

("C/Wm-2) 

circuit  and  the  ambient  (°C) 

Rtb 

absolute  thermal  resistance  between  the  solar  cells  in 

At mpp 

temperature  difference  between  the  cell  at  maximum 

a  module  and  the  ambient  (°C/W) 

power  point  and  the  ambient  (°C) 

Tair 

ambient  temperature  (°C) 

fjmeas 

averaged  electrical  efficiency  of  an  HCPV  module 

Tc 

cell  temperature  (°C) 

while  measuring  an  l-V  curve  (dimensionless) 

Tc(K) 

absolute  cell  temperature  (K) 

Rmpp 

electrical  efficiency  of  an  HCPV  module  working  at 

W 

reference  absolute  cell  temperature  (K) 

maximum  power  point  (dimensionless) 

Tc.meas 

cell  temperature  measured  at  open-circuit  (°C) 

lopt 

optical  efficiency  (dimensionless) 

Tc,mpp 

cell  temperature  working  at  maximum  power  point 

A 

thermal  conductivity  of  a  material  (W  m_1/°C) 

(°C) 

P 

internal  thermal  resistance  of  an  HCPV  module 

Th-s 

module  heat-sink  temperature  (°C) 

(°C/W  m~2) 

Th-s,o 

heat-sink  temperature  at  shutter  initiation  (DC) 

Qi,  a2,  b1 

,  b2  parameters  of  the  method  of  Peharz  (V,  V,  V/°C, 

Vac 

open-circuit  voltage  (V) 

V/“C) 

V  * 
v  oc 

open-circuit  voltage  at  reference  conditions  (V) 

Ai,  A2,  A3 parameters  of  the  atmospheric  method  for  heat-sink 

V oc, max 

maximum  open-circuit  voltage  registered  during 

temperature  calculation  of  PVUSA  (dimensionless, 

a  shuttering  procedure  (V) 

°C/W  m~2,  °C/m  s_1) 

ws 

wind  speed  (m/s) 

Amod 

module  aperture  area  (m2) 

WSo 

parameter  of  the  method  of  Hornung  (m/s) 

Bj,  B 2 

parameters  of  the  atmospheric  method  for  heat-sink 

Xc 

concentration  ratio  (dimensionless) 

temperature  calculation  of  Sandia  (dimensionless, 

1/m  s_1) 

or  mechanics,  that  allows  displacing  the  manufacturing  process  from 
the  pure  semiconductor  industry.  This  is  beneficial  for  a  fast  growth  of 
this  technology  and  offers  opportunities  for  costs  reduction  [28],  HCPV 
systems  mainly  combine  the  following  elements. 


2.1.  Solar  cells 

The  photovoltaic  cells  used  in  these  systems  are  designed  to 
extract  more  current  per  unit  of  area  than  conventional  cells.  Solar 


cell  open-circuit  voltage  and  fill  factor  grow  with  the  incident 
luminous  flux  and,  as  a  consequence,  conversion  efficiency 
increases  the  more  sunlight  is  concentrated  (until  a  limit  given 
by  the  Ohmic  losses).  Silicon-based  solar  cells  can  be  used  in 
HCPV,  but  III— V  semiconductors-based  cells  are  more  suitable  for 
concentration  because  of  the  lower  base  electrical  resistivity. 
Nowadays,  advanced  III— V  multi-junction  solar  cells  have  demon¬ 
strated  the  highest  efficiencies  because  of  a  broad  use  of  the 
sunlight  spectrum  and  many  HCPV  systems  are  based  on  these 
kinds  of  cells  [29], 
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increases. 
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Fresnel  lens 
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Fig.  3.  Structure  of  a  solar  receiver  integrated  with  Fresnel  lens-based  primary 
concentrator. 


Classification  of  methods  for 
cell  temperature  calculation 
in  HCPV  modules 


Methods  based  on  module 
heat-sink  temperature 

Methods  based  on  direct 
measurements  on  the  module 

Methods  based  on  electrical 
parameters 


Methods  based  on 
atmospheric  parameters 


Fig.  2.  Classification  of  the  reviewed  methods  for  the  calculation  of  cell  tempera¬ 
ture  in  HCPV  modules. 


2.2.  Optics 

Lenses  and/or  mirrors  are  used  to  increase  the  luminous  flux  on 
the  solar  cells.  Primary  optics  and  secondary  optics  can  be 
distinguished.  Primary  optics  is  intended  for  obtaining  a  high 
concentration  level.  Parabolic  mirrors  or  Fresnel  lenses  are  com¬ 
monly  used  [30],  Secondary  optics  is  optional.  It  receives  the  light 
from  the  primary  concentrator  and  focuses  it  on  the  solar  cell  with 
the  aim  of  light  homogenization  and  angular  acceptance  improve¬ 
ment  [31,32], 

Several  magnitudes  related  to  the  optics  in  HCPV  are 

The  geometric  concentration  ( Cg ).  It  is  the  quotient  between  the 
primary  concentrator  area  and  the  solar  cell  area. 

The  optical  efficiency  ( rjopt ).  It  is  the  quotient  between  the 
irradiation  actually  incident  on  the  solar  cell  and  the  collected 
irradiation.  As  HCPV  systems  only  exploit  the  direct  component 
of  sunlight,  the  collected  irradiation  is  the  product  of  the  direct 
normal  irradiance  (DM)  by  the  primary  concentrator  area.  r/op[ 
accounts  for  the  optical  losses  due  to  non-ideal  properties  of 
reflection,  transmission  and  absorption  of  the  optical  elements. 
The  concentration  ratio  (Xc).  It  is  the  quotient  between  the  mean 
irradiance  actually  incident  on  the  solar  cell  and  the  one-sun 
irradiance  (1000  W/m2). 


2.3.  Solar  receivers 

A  solar  receiver  is  an  assembly  of  one  or  several  solar  cells  with 
the  secondary  optics  and  the  mechanisms  for  current  extraction 
and  heat  dissipation.  It  can  incorporate  other  elements  such  as  a 
bypass  diode  in  order  to  prevent  cell  overheating  in  shading 
conditions  (in  two-axis  trackers  plants,  modules  on  a  tracker  are 
shaded  by  the  neighboring  trackers  at  certain  hours  of  the  day) 
[33,34],  Fig.  3  shows  the  structure  of  a  typical  solar  receiver 
integrated  with  a  Fresnel  lens-based  primary  optics.  In  this  case, 


passive  cooling  is  achieved  by  means  of  a  finned  heat  exchanger. 
Some  of  the  methods  reviewed  in  this  paper  require  measuring 
the  heat-sink  temperature  (Th_s)  by  placing  a  temperature  sensor 
at  the  root  of  the  heat  exchanger  behind  the  solar  cell. 


2.4.  HCPV  modules 

An  HCPV  module  is  the  smallest,  complete,  environmentally 
protected  assembly  of  solar  receivers  and  optical  devices  that  is 
able  to  transform  an  input  of  unconcentrated  solar  radiation.  Fig.  4 
shows  the  main  components  and  the  structure  of  a  typical  HCPV 
module  based  on  Fresnel  lens  parquets. 


2.5.  Tracking  system 

HCPV  modules  are  mounted  on  a  mobile  structure  in  order  to 
be  always  oriented  to  collect  the  solar  rays.  Some  linear-focus 
systems  can  operate  on  single-axis  solar  trackers.  However,  two- 
axis  trackers  are  the  most  used  and  are  required  for  point-focus 
systems  [35,36].  These  trackers  are  more  complex  from  a  mechan¬ 
ical  point  of  view  but  offer  the  advantage  of  higher  light  concen¬ 
tration.  Solar  trackers  are  composed  of  a  tracking  structure 
(structural  elements  and  moving  mechanisms)  and  the  tracking 
control  (electronic  equipment  that  governs  the  motion).  A  set  of 
modules  mounted  on  the  same  tracker  is  known  as  a  concentrator 
array,  while  a  set  of  concentrator  arrays  connected  to  the  same 
power  output  constitute  an  HCPV  field  [28], 


3.  Methods  based  on  module  heat-sink  temperature 

In  this  section,  three  methods  for  cell  temperature  calculation 
are  summarized  that  require  measuring  the  module  heat-sink 
temperature  (Th_s).  The  heat-sink  temperature  is  typically  mea¬ 
sured  at  the  root  of  a  finned  heat  exchanger  (heat  sink)  on  the  back 
of  the  module.  This  temperature  is  distinctly  different  than  the  cell 
temperature  in  the  case  of  HCPV  modules  and  mainly  depends  on 
direct  normal  irradiance,  position  of  the  sensor,  module  design, 
ambient  temperature,  wind  speed  and  direction  and  module 
electrical  efficiency  [26,37],  The  direct  measurement  of  the  heat¬ 
sink  temperature  avoids  using  thermal  models  that  try  to  estimate 
this  temperature  from  atmospheric  parameters  but  present 
uncertainty. 
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MODULE 


RECEIVER 


Fig.  4.  Structure  of  an  HCPV  module  based  on  point-focus  Fresnel  lenses. 
Source:  IEC  62108  [1  ]. 


3.1.  Method  of  Sandia's  photovoltaic  array  performance  model 
(2004) 

The  method  presented  by  King  et  al.  [38]  tries  to  express  the 
cell  temperature  as  a  function  of  the  heat-sink  temperature  and 
the  DM.  For  this  purpose,  it  assumes  one-dimensional  thermal 
heat  conduction  through  the  module  materials  behind  the  cell.  It  is 
based  on  a  simple  equation  as  follows: 

Tc  =  Th_s  +  (DNI/DNI*)AT  (1) 

AT  being  the  temperature  difference  between  the  cell  and  the 
heat-sink  at  the  defined  reference  direct  normal  irradiance  DM*. 
Obtaining  the  AT  parameter  is  not  easy  because  it  requires 
accessing  inside  the  module  in  order  to  measure  the  cell  tempera¬ 
ture  at  the  reference  irradiance. 

3.2.  Method  of  ISFOC/IES-UPM  (2008) 

The  method  developed  by  the  Institute  de  Sistemas  Fotovoltaicos 
de  Concentration  (ISFOC)  in  collaboration  with  the  Institute  de  Energia 
Solar  of  the  Universidad  Politecnica  de  Madrid  (IES-UPM)  [39]  uses  a 
similar  assumption  than  the  previous  one,  i.e.  the  relation  between  the 
cell  and  the  heat-sink  temperature  is  expressed  as  a  function  of  DM: 

Tc  =  Ths+pDNI  (2) 

The  main  difference  with  respect  to  the  previous  method  is  that 
the  equation  depends  on  a  physical  parameter,  the  so  called 
“internal  thermal  resistance”  ( p ,  in  °C/W  m~2),  which  is  propor¬ 
tional  to  the  thermal  resistance,  R,  between  the  cell  and  the  heat¬ 
sink.  The  p  parameter  can  be  theoretically  calculated  if  the 
detailed  structure  of  the  materials  behind  the  cell  is  known,  as 
well  as  the  geometric  concentration  ratio,  Cg,  and  optical  effi¬ 
ciency,  t/0pt,  of  the  module.  Proposed  expressions  for  calculating  p 
are  [40] 

K=ia,A)  (3) 

i 

P  =  >1optcgR  (4) 

where  index  i  represents  each  layer  of  material  behind  the  cell, 
U  and  2,-  being  the  thickness  and  thermal  conductivity  of  the 
material  respectively.  The  theoretical  calculation  of  p  avoids  the 
need  of  accessing  inside  the  module;  however,  a  detailed 


knowledge  of  the  module  design  is  required  to  perform  this 
calculation  and  this  information  is  not  always  available. 

3.3.  Method  of  Muller  et  al.  (2011 ) 

The  method  presented  by  Muller  et  al.  [41]  from  the  National 
Renewable  Energy  Laboratory  (NREL)  requires  measuring  both  the 
heat-sink  temperature  and  the  open-circuit  voltage.  It  has  been 
included  into  this  category  of  methods  because  it  requires  the 
heat-sink  temperature,  although  it  also  could  have  been  classified 
as  a  method  based  on  electrical  parameters.  The  method  tries  to 
overcome  the  limitations  of  common  methods  based  on  the  open- 
circuit  voltage,  which  are  very  sensitive  to  certain  parameters  not 
always  easily  obtainable,  such  as  the  effective  diode  ideality  factor 
of  the  module. 

One  of  the  difficulties  of  this  method  is  that  it  requires  a  system 
that  shutters  sunlight  to  the  module  while  measuring  the  tran¬ 
sients  in  open-circuit  voltage.  In  order  to  get  accurate  results,  a  fast 
shutter  must  be  available  because  the  procedure  is  sensitive  to  the 
aperture  time  as  well  as  a  fast  multimeter  in  order  to  be  able  to 
register  the  transients  in  open-circuit  voltage.  Also,  the  module 
must  be  covered  for  a  significant  period  of  time  before  all 
shuttering  events  in  order  to  equilibrate  the  cell  and  heat-sink 
temperatures  at  the  beginning  of  the  procedure. 

When  the  procedure  is  launched,  the  open-circuit  voltage 
quickly  grows  until  reaching  a  maximum  ( Vocmax ).  Afterwards,  it 
gradually  decreases  as  the  cell  temperature  increases  until  equili¬ 
brium  is  reached.  The  heat-sink  temperature  grows  with  time 
towards  the  equilibrium  state  and  so  does  the  cell  temperature. 

The  following  equation  is  proposed  to  model  the  behavior  of 
the  cell  temperature  during  transient  at  an  instant  t: 

Tc(t)  =  T/i  _  so  —  [V  oc  max-  Voc{t)\/f)Voc  (5) 

T),_Sio  being  the  measured  heat-sink  temperature  at  shutter  initia¬ 
tion  and  pVoc  the  temperature  coefficient  of  Voc  for  the  module.  The 
cell  operating  temperature  can  be  obtained  by  evaluating  this 
expression  at  an  instant  t  in  which  the  measured  open-circuit 
voltage  can  be  considered  stable. 

Authors  of  the  method  compared  some  results  against  manu¬ 
facturers'  suggestions  and  against  a  conventional  method  based  on 
the  open-circuit  voltage.  Results  were  in  agreement  for  fast 
shutters.  However,  an  experimental  validation  based  on  compar¬ 
ing  against  measured  cell  temperatures  has  not  yet  been  reported. 
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4.  Methods  based  on  electrical  parameters 

In  this  section,  the  methods  based  on  electrical  parameters  are 
summarized.  These  methods  are  based  on  the  fact  that  the  cell 
electrical  parameters,  specially  the  open-circuit  voltage,  vary  with 
temperature.  This  allows  estimating  the  cell  temperature  from 
the  measurement  of  electrical  parameters  by  using  different 
procedures. 

4.1.  Method  of  the  standard  IEC  60904-5  (2011 ) 

The  method  of  the  standard  IEC  60904-5  [42]  is  known  as  the 
open-circuit  voltage  method.  It  requires  measuring  the  module 
open-circuit  voltage  together  with  the  incident  DM.  The  measured 
open-circuit  voltage  ( Voc )  can  be  related  with  the  open-circuit 
voltage  at  reference  conditions  ( Voc *)  as  a  function  of  absolute  cell 
temperature  (T^k))  and  DM 


Voc  =  V„c*  +  Ns(nkTc(K)/q)\n(DNI/DN!*)  +fJVoc(Tc(K)  -  TC(K)*)  (6) 

Ns  being  the  number  of  cells  in  series  in  the  module,  k  the 
Boltzmann  constant  and  q  the  electron  charge.  The  defined 
reference  conditions  of  irradiance  and  absolute  cell  temperature 
are  denoted  by  DM*  and  Tc(Kl*  respectively.  The  expression  also 
depends  on  two  important  parameters:  the  temperature  coeffi¬ 
cient  of  Voc  for  the  module  (/W)  and  the  effective  diode  ideality 
factor  of  the  module  (n). 

From  this  equation,  the  absolute  cell  temperature  can  be 
expressed  as 

Tc(K)  =  (Voc-V0c*+/3VocTciK*)/[Ns(nk/q)ln(DNI/DNI*)+/3Voc]  (7) 

Expression  that  allows  calculating  the  cell  temperature  as  a 
function  of  Voc  and  DM. 

While  the  f)Voc  parameter  is  sometimes  provided  by  the 
manufacturer,  the  n  parameter  is  not  easily  obtainable.  If  these 
values  are  taken  from  the  literature,  there  is  the  problem  of 
uncertainty  because  the  method  is  very  sensitive  to  the  chosen 
values. 

A  way  of  overcoming  these  problems  is  to  carry  out  an 
outdoor  experimental  campaign  in  order  to  register  measure¬ 
ments  of  cell  temperature,  open-circuit  voltage  and  DM  over  a 
significant  period  of  time.  For  this  purpose,  a  module  with  a  cell 
temperature  sensor  must  be  available.  From  the  collected  data,  a 
multi-linear  regression  analysis  can  be  done  taking  into  account 
Eq.  (6)  and  a  better  estimate  of  the  pVoc  and  n  parameters  can  be 
achieved  [40], 

An  alternative  to  this  method  has  been  proposed  by  King  et  al.  [38] 
that  avoids  the  use  of  a  pyrheliometer.  It  consists  of  approximating  the 
DM/DM*  quotient  by  the  relation  between  the  measured  short-circuit 
current  (/5C)  and  the  short-circuit  current  at  reference  conditions  (/sc*). 
This  alternative  is  known  as  the  “Voc-4c  method"  [41,43]  and  thus,  is 
based  on  the  following  equation: 

Tc(K)  =  (Voc-Voc*+/JVocTc(K*)/[Ns(nk/q)\n(Isc/Isc*)+PVoc]  (8) 


4.2.  Method  ofju  et  al.  (2013) 

The  method  proposed  by  Ju  et  al.  [43  can  be  regarded  as  an 
improvement  of  the  Voc-Isc  method,  allowing  the  variation  of  the 
temperature  coefficient  of  Voc  (fivoc)  with  the  irradiance  level. 
Several  authors  have  shown  that  the  temperature  coefficient  of  Voc 
for  a  triple  junction  III— V  solar  cell  ( /W.ce/f )  can  be  expressed  as  a 
logarithmic  function  of  the  concentration  ratio,  Xc,  in  suns  [7,44] 


Pvoc ,cei/(Xc=  1 )  being  the  temperature  coefficient  of  Voc  for  the  cell 
at  one  sun  and  nceU  the  diode  ideality  factor  of  the  cell. 

Based  on  this  equation,  it  can  be  assumed  that  the  temperature 
coefficient  for  an  HCPV  module  could  be  approximated  by 

fivoc  fivocO  +Ns{nk/q)\n(DNI/DNl*) «  fiVocQ +Ns(nk/q)\n{Isc/Isc*) 

(10) 


fivoco  being  the  temperature  coefficient  of  Voc  for  the  module  at  the 
defined  reference  irradiance  DM*. 

Taking  this  into  account,  the  method  is  derived  from  Eq.  (6)  by 
neglecting  the  variations  of  the  thermal  voltage  term  over  the 
range  of  operation  temperatures.  This  equation  can  be  rewritten  as 

Voc  =  Voc*+Ns(nkTc(K)*/q)ln(Isc/IsC*)+/3Voc(TC(K)  -  W)  (11) 


By  replacing  Eq.  (10)  into  Eq.  (11),  an  expression  for  calculating  the 
absolute  cell  temperature  is  achieved 


j.  _  ..  *  Voc-  V oc  —Ns(nkTC(K)  /q)ln(Isc/Isc  ) 

ciK)~  aK)  flVoco+Ns(nk/q)ln(Isc/Isc*) 


(12) 


The  method  depends  on  the  same  parameters  than  the  Voc-Isc 
method  so  that  it  can  be  used  as  an  alternative  that  theoretically 
could  improve  the  estimation.  However,  it  is  important  to  remark 
that  the  method  was  not  validated  by  the  authors  for  HCPV 
modules  operating  outdoors,  so  that  the  usefulness  of  the  method 
for  the  case  of  HCPV  modules  is  still  unknown.  Authors  evaluated 
the  method  for  a  single  cell  measured  in  a  solar  simulator  under 
different  concentration  levels  ranging  from  1  to  1000  suns.  From 
these  measurements,  the  method  presented  better  results  than  the 
Voc-Isc  method  especially  at  high  concentrations. 


4.3.  Method  of  Peharz  et  al.  (2011 ) 

The  method  presented  by  Peharz  et  al.  [18]  from  the  Fraunhofer 
Institut  ftir  Solare  Energiessysteme  is  based  on  an  indoor  char¬ 
acterization  procedure  applied  to  several  HCPV  modules.  Different 
measurements  carried  out  in  a  sun  simulator  allowed  the  authors 
deriving  a  functional  dependence  of  the  cell  temperature  from  the 
module  electrical  parameters  of  open-circuit  voltage  and  short- 
circuit  current. 

The  sun  simulator  used  for  the  study  [45]  allows  controlling  the 
light  intensity  incident  on  the  module  by  means  of  a  flash  bulb. 
The  module  temperature  can  also  be  controlled  by  irradiating  the 
module  backside  by  several  infrared  light  bulbs.  This  temperature 
can  be  varied  by  changing  the  distance  between  the  bulbs  and  the 
module.  After  about  30  min  of  infrared  irradiation,  the  back 
temperature  of  the  module  is  stable  and  can  be  measured  with  a 
thermocouple  placed  on  uniformly  spaced  points,  what  allows 
determining  the  average  temperature  across  the  backside  of  the 
module.  A  theoretical  analysis  based  on  Fourier's  heat  conduction 
equation  showed  that  this  average  temperature  is  almost  the  same 
than  the  cell  temperature  at  the  described  thermal  equilibrium 
condition  (with  a  difference  less  than  0.1  °C).  Thus,  the  cell 
temperature  can  be  controlled  in  the  simulator  and  can  be 
determined  from  the  measurements  of  the  thermocouple. 

With  this  experimental  set-up,  the  Voc  and  Isc  of  the  modules 
under  study  were  systematically  measured  at  different  conditions 
of  irradiance  and  temperature.  It  was  found  that  the  open-circuit 
voltage  fitted  the  following  function: 

V0c  =  [Qi  lti(/sc) -f  U2] 3“ [hi  ln(/sc) +  b2]fc  (33) 

and  from  this  equation,  the  cell  temperature  can  be  expressed  as 

T  _  Voc-Qi  ln(fsc)-a2 

c~  bj  ln(Isc)+b2 


Pvoc.ceii  —  Pvoc.ceiiWc  —  1 ) +(ncenk/q)\n(Xc) 


0) 


(14) 
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The  main  difficulty  of  the  method  is  that  a  complex  experi¬ 
mental  set-up  is  required  for  obtaining  the  alt  a2,  tq  and  b2 
parameters  for  a  particular  HCPV  module.  However,  in  the  present 
research,  we  propose  an  alternative  to  this  procedure  based  on 
outdoor  measurements.  It  consists  on  monitoring  the  HCPV 
module  during  an  experimental  campaign  of  several  months  in 
order  to  register  measurements  of  cell  temperature,  open-circuit 
voltage  and  short-circuit  current  over  a  representative  range  of 
values.  From  the  measurements,  a  multi-linear  regression  analysis 
can  be  done  taking  into  account  Eq.  (13)  that  allows  the  determi¬ 
nation  of  the  four  parameters.  Results  of  this  outdoor  procedure  is 
shown  in  Section  6. 


where  Pmpp(Tc,mpp)  is  the  power  of  the  MPP  module  operating  at 
Tc,mpp  temperature  and  PmeasAvc(Tc.meas )  is  the  averaged  power 
harvested  by  the  measured  module  during  testing  (while  perform¬ 
ing  the  voltage  sweep  needed  to  trace  the  I-V  curve)  at  TCimeas 
temperature.  The  averaged  power  harvested  by  the  measured 
module  can  be  calculated  from  the  registered  P-V  curve: 

rVOC 

P  meas/iVG  =  (l/v0c)  /  P(V)dV  (22) 

Jo 

Thus,  from  Eq.  (19)  the  following  can  be  written: 

AT  mpp  _  A  Tmeas  (23) 

1  ~  P MPp(TcMPp) / Tabs  1  ~  PmeasAVG(T c.meas) / P abs 


4.4.  Method  of  Fernandez  et  al.  (2013) 

The  method  proposed  by  Fernandez  et  al.  [46]  uses  linear 
coefficients  in  order  to  express  the  dependence  of  open-circuit 
voltage  (Voc)  on  DN1  and  cell  temperature  (Tc) 

V  oc  ==  C\DNI  -f-  c2T  c  -f  Cj  (15) 

This  way,  once  the  coefficients  have  been  determined,  the  cell 
temperature  is  calculated  by 

Tc  =  (V0C-CiDNI-c3)/c2  (16) 

The  method  is  simple  but  requires  an  outdoor  experimental 
campaign  on  a  monitored  module  in  order  to  obtain  the 
coefficients. 


4.5.  Method  of  Yandt  et  al.  (2012 ) 

The  above  methods  are  intended  for  calculating  the  cell 
temperature  in  an  HCPV  module  at  open-circuit  condition. 
An  interesting  aspect  in  HCPV  is  to  investigate  the  cell  tempera¬ 
ture  of  the  module  working  at  maximum  power  point  (MPP).  At 
MPP,  the  module  is  harvesting  electricity  and  this  energy  is  not 
transformed  to  heat.  As  a  result,  the  cell  temperature  at  MPP 
( Tc,mpp )  will  be  lower  than  the  cell  temperature  as  measured  at 
open-circuit  ( Tcmeas ).  Knowing  the  cell  temperature  at  MPP  can  be 
interesting  in  certain  studies,  for  instance  in  the  energy  harvesting 
calculation  of  an  HCPV  system. 

Yandt  et  al.  [47]  proposed  a  method  for  calculating  the  cell 
temperature  at  MPP.  It  is  based  on  correcting  Tcmeas  in  order  to 
estimate  TcMPP.  The  method  requires  measuring  the  whole  I-V 
curve  of  a  monitored  module  and  uses  the  atmospheric  conditions 
of  ambient  temperature  (Tair)  and  DM. 

The  absolute  thermal  resistance  ( Rth )  between  the  solar  cell  and 
the  ambient  is  used  to  derive  the  method.  This  magnitude  can  be 
defined  as 


AT, 


=  AT 


meas 


1  —  1]  mpp(Tc.MPp) 

1  ~0meas(T  c.meas ) 


(24) 


>1mpp'  Omeas  being  the  efficiencies  related  to  the  luminous  power 
absorbed  by  the  cells  for  the  cases  of  the  MPP  module  and  the 
measured  module  respectively. 

In  Eq.  (24),  Ompp(Tc.mpp)  is  unknown.  However,  a  first  approach 
to  this  value  is  OMPtfTc.meas ),  which  can  be  extracted  from  the 
maximum  power  of  the  measured  I-V  curve.  Taking  this  into 
account,  the  authors  propose  an  iterative  method  for  calculating 
A Tmpp.  It  is  based  on  the  following  sequence: 


AT, 


MPP(n) 


=  A  Tn 


1  0 MPP^T c,MPP{n  -  1)) 

1  ~  0 meas(T c.meas) 


(25) 


where  TCiMpp(o)=Tc,meas  and  ATmpp=  ATmpKoo)- 

One  of  the  difficulties  of  the  method  is  that  the  luminous 
power  absorbed  by  the  cells  is  not  easily  obtainable.  A  first 
approach  to  this  power  could  be 


Pabs  ~~  OopfiNIAmod 


(26) 


Oopt  being  the  module  optical  efficiency  and  Amod  the  module 
aperture  area.  However,  not  all  of  this  quantity  is  actually  absorbed 
by  the  cell  because  some  light  is  reflected  and  some  light  is 
transmitted  [48-50],  Also,  as  it  is  known,  rjopt  depends  on  factors 
such  as  incident  spectrum  and  lens  temperature,  so  that  establish¬ 
ing  a  constant  value  for  the  optical  efficiency  would  be  only  an 
approximation  [51-53], 

It  is  remarkable  that  the  authors  of  the  method  have  not  yet 
published  experimental  validation,  i.e.  results  of  the  method 
compared  against  measured  cell  temperatures  at  MPP  have  not 
yet  been  reported. 


5.  Methods  based  on  atmospheric  parameters 


Pth  —  (T c  ~  T a,,- / Q)  (12) 

Q  being  the  rate  of  heat  transfer  through  conduction  and  convec¬ 
tion.  This  definition  can  be  applied  to  both  the  MPP  module  and 
the  measured  module 

Pth  —  (Tc.MPP  ~  Tair) / Q.MPP  —  (T  c.meas  —  T  a,r)  /  Qmeas  (18) 

Pth  =  ATmpp/Qmpp  =  AT  meas  /Q.  meas  (19) 

The  Qmpp  and  Q,,1(.as  rates  can  be  expressed  as  the  difference 
between  the  luminous  power  absorbed  by  the  cells  in  the  module 
(Pabs)  and  the  electrical  power  harvested  by  the  module 

0.MPP  =  Pabs  —  PmPp(Tc.MPp)  (20) 

(21) 


The  methods  reviewed  above  require  some  direct  measure¬ 
ments  on  a  module  so  that  they  are  adequate  during  field 
performance  studies.  However,  when  the  set-up  needed  for  the 
measurements  is  not  available  or  when  dealing  with  studies  based 
on  typical  meteorological  data,  that  practice  is  not  possible. 
In  these  cases,  methods  only  based  on  atmospheric  parameters 
must  be  used  for  estimating  the  cell  temperature. 

An  important  task  for  the  analysis  of  an  HCPV  installation  is  to 
estimate  the  energy  harvesting.  This  is  done  from  values  of  several 
atmospheric  parameters,  such  as  DM,  ambient  temperature  or 
wind  speed,  among  others.  As  the  cell  temperature  is  an  important 
factor  for  energy  harvesting  calculations,  it  is  very  interesting  to 
have  methods  that  exclusively  estimate  this  temperature  from 
atmospheric  parameters.  These  procedures  present  uncertainties 
related  to  the  sources  of  the  meteorological  data  and  related  to  the 
model  itself  but  have  proven  to  be  a  very  useful  tool  [54], 


Q.meas — Pabs  PmeasAVG(Tc.meas) 
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5.1.  Method  of  Almonacid  et  al.  (2012) 

The  method  presented  by  Almonacid  et  al.  [13]  uses  linear 
coefficients  in  order  to  estimate  the  cell  temperature  in  an  HCPV 
module  from  ambient  temperature  (Tair),  DM  and  wind  speed 
(WS).  It  is  based  on  the  following  relation: 

Tc  =  Tair  +  d1DNI+d2WS  (27) 

where  dr  and  d2  parameters  must  be  obtained  by  means  of  multi¬ 
linear  regression  analysis  of  monitored  data  for  the  analyzed  HCPV 
module. 

5.2.  Method  of  Hornung  et  al.  (2012) 

A  slightly  more  sophisticated  method  was  proposed  by  Hor¬ 
nung  et  al.  [55]  based  on  the  same  atmospheric  parameters  than 
the  previous  one.  In  this  case,  the  cell  temperature  is  estimated  by 

Tc  =  Tair  +  m[exp(-0.5WS/WSo)  +  c]DM  (28) 

The  method  assumes  the  temperature  gap  between  the  cell  and 
the  ambient  to  be  proportional  to  the  DM,  with  an  exponential 
correction  based  on  wind  speed.  A  multi-linear  regression  analysis 
can  be  done  to  obtain  the  model  coefficients  from  monitored  data. 

5.3.  Method  of  Fernandez  et  al.  based  on  ANN  (2013) 

The  above  two  methods  use  simple  algebraic  functions  for 
characterizing  the  relation  between  cell  temperature  and  atmo¬ 
spheric  parameters.  Actually,  this  relation  is  more  complex. 
Because  of  this,  Fernandez  et  al.  [40]  tried  to  characterize  this 
relation  by  using  an  artificial  neural  network  (ANN)  in  order  to 
find  a  more  accurate  model.  ANNs  have  proven  to  be  very  useful 
for  solving  complex  problems  and  have  been  recently  applied  in 
the  HCPV  field  [56,57], 


The  same  atmospheric  parameters  than  the  previous  methods 
are  proposed  to  characterize  the  cell  temperature,  i.e.  the  objective 
is  to  find  the  function 

Tc=f(Tair.  DN1.WS)  (29) 

To  find  this  function,  a  feed-forward  ANN  composed  of  three 
layers  was  implemented:  a  three  nodes  layer  ( Tair ,  DM  and  WS)  as 
input  layer,  a  five  nodes  layer  as  hidden  layer  and  a  single  node 
layer  (Tc)  as  output  layer.  The  number  of  nodes  of  the  hidden  layer 
was  empirically  determined.  This  architecture  was  trained  with 
the  Levenberg-Marquardt  back-propagation  algorithm,  which 
allows  determining  the  coefficients  of  the  ANN  (weights  and 
biases  of  the  hidden  and  output  layers).  In  order  to  train,  validate 
and  test  the  ANN,  a  set  of  outdoor  measurements  including  the 
cell  temperature  and  the  atmospheric  parameters  must  be  avail¬ 
able  for  a  wide  range  of  operating  conditions. 

The  main  difficulty  of  the  method  is  that  an  advanced  knowl¬ 
edge  on  ANNs  is  required  to  train  the  network.  However,  once  the 
coefficients  of  the  network  are  determined,  the  method  becomes 
easy-to-use. 

After  completing  the  review  of  existing  methods  for  the 
calculation  of  cell  temperature  in  an  HCVP  module,  it  is  important 
to  remark  that  some  efforts  have  been  made  in  order  to  be  able  to 
estimate  the  module  heat-sink  temperature  from  atmospheric 
parameters.  These  methods  can  also  be  useful  when  the  direct 
measurement  of  this  temperature  on  the  module  is  not  possible. 
The  application  of  these  methods  allows  converting  the  methods 
based  on  module  heat-sink  temperature  (Section  3)  to  methods 
only  based  on  atmospheric  parameters.  Two  proposed  expressions 
are: 

-  From  photovoltaics  for  utility  scale  applications  (PVUSA)  [58] 

Th_s=ArTair+A2DNI+A3WS  (30) 


Table  1 

Summary  of  characteristics  of  the  methods  for  calculating  the  cell  temperature  in  an  HCPV  module:  required  inputs,  advantages  and  technical  difficulties. 


Methods  based  on 

Methods  based  on  electrical  parameters 

Methods  based  on  atmospheric 

module  heat-sink 
temperature 

parameters 

Sandia 

ISFOC  / 

Muller 

IEC 

Ju  Peharz 

Fernandez 

Yandt 

Almonacid 

Hornung 

Fernandez 

IES- 

60904-5 

ANN 

UPM 

Required  inputs 

Module  heat-sink  temperature 

V 

V 

V 

Direct  normal  irradiance 

V 

V 

V 

V 

V 

V 

V 

V 

Open-circuit  voltage 

V 

V 

V  V 

V 

Short-circuit  current 

Cell  temperature  at  open-circuit 

V  V 

V 

Whole  I-V  curve 

V 

Ambient  temperature 

V 

V 

V 

V 

Wind  speed 

V 

V 

V 

Advantages 

Outdoor  validation  has  been  reported 

V 

V 

V 

V 

V 

V 

V 

V 

Parameters  can  be  obtained  from  manufacturers 

V 

V 

V 

V 

V 

or  from  the  literature 

Does  not  require  direct  measurements  on  the  module 
Intended  for  calculating  cell  temperature  at  MPP 

V 

V 

V 

V 

Technical  difficulties 

Requires  module  with  cell  temperature  sensor  to  be 

V 

V 

V 

V 

V 

measured  outdoors  for  determining  the  parameters 
Requires  a  sun  simulator  for  determining  the  parameters 
Requires  a  system  that  shutters  sunlight 

Requires  the  whole  I-V  curve  of  the  module 

Requires  advanced  knowledge  on  ANN 

V 

Va 

V 

V 

Alternatively,  parameters  can  be  determined  from  outdoor  measurements  on  a  module  with  cell  temperature  sensor. 
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-  From  Sandia  National  Laboratories  [38] 

Th-s  =  rm-r  +  exp(Bi  +B2WS)DNI.  (31) 


6.  Comparative  analysis 

Table  1  summarizes  the  required  inputs,  advantages  and 
technical  difficulties  of  the  reviewed  methods  for  calculating  the 
cell  temperature  in  an  HCPV  module.  As  can  be  seen,  methods 
based  on  direct  measurements  on  the  module  are  preferable  for 
field  performance  analysis  while  methods  based  on  atmospheric 
parameters  are  suitable  when  direct  measurements  on  the  module 
are  not  possible. 

Among  the  methods  based  on  direct  measurements,  some  of 
them  have  the  advantage  that  the  parameters  can  be  obtained  from 
manufacturer's  data  or  from  the  scientific  literature  (ISOFC/IES-UPM, 
Muller,  IEC  60904-5,  Ju  and  Yandt).  However,  as  has  been  commen¬ 
ted,  accurate  values  of  the  parameters  are  required  because  these 
methods  are  very  sensitive  to  the  chosen  values.  The  rest  of  the 
methods  based  on  direct  measurements  have  the  difficulty  that 
parameters  must  be  fitted  from  an  outdoor  experimental  campaign 
or  from  indoor  measurements  in  a  sun  simulator. 

Methods  based  on  atmospheric  parameters  always  require  an 
outdoor  experimental  campaign  in  order  to  obtain  the  parameters 
and  use  values  of  direct  normal  irradiance,  ambient  temperature 
and  wind  speed  as  inputs. 

The  table  remarks  some  technical  difficulties  related  to  the 
application  of  the  methods.  As  has  been  commented,  some  methods 
require  an  outdoor  experimental  campaign  on  a  module  with  a  cell 
temperature  sensor  in  order  to  determine  the  parameters.  This  kind  of 
device  is  not  often  available  and  the  procedure  takes  time  to  carry  out 
the  measurements.  Other  technical  difficulties  are  the  need  of  a  sun 
simulator  (Peharz),  a  system  that  shutters  sunlight  (Muller),  the  need 
of  the  whole  I-V  curve  of  the  module  (Yandt)  or  the  need  of  advanced 
knowledge  in  ANN  (Fernandez  ANN). 

As  can  be  seen  in  the  table,  only  one  method  is  intended  for 
calculating  the  cell  temperature  of  a  module  operating  at  max¬ 
imum  power  point  (Yandt).  Authors  of  this  method  predicted  a 
difference  of  up  to  14  C  between  the  cell  temperature  at  open- 
circuit  and  the  cell  temperature  at  MPP  [47],  Thus,  investigating 
the  cell  temperature  at  MPP  is  interesting  in  HCPV  and  further 
research  in  this  field  is  expected.  In  addition,  it  is  important  to 
remark  that  there  are  three  methods  whose  outdoor  validation 
against  measured  cell  temperature  has  not  yet  been  reported 
(Muller,  Ju  and  Yandt). 

A  quantitative  comparison  of  the  methods  has  also  been  carried 
out.  For  this  purpose,  an  HCPV  module  equipped  with  cell 
temperature  sensor  and  heat-sink  temperature  sensor  was  mea¬ 
sured  for  two  years,  from  January  2011  to  December  2012,  at  Jaen, 
South  of  Spain.  The  temperatures  were  recorded  every  5  min 
together  with  the  main  atmospheric  parameters  (direct  normal 
irradiance,  ambient  temperature  and  wind  speed).  Also,  an  I-V 
tracer  simultaneously  measured  the  module  I-V  curve,  from  which 
the  electrical  parameters  open-circuit  voltage  and  short-circuit 
current  can  be  extracted.  The  main  characteristics  of  the  measured 
module  are  depicted  in  Table  2. 


In  order  to  carry  out  the  experimental  campaign,  the  module 
was  mounted  on  a  two-axis  solar  tracker  located  at  the  Centro  de 
Estudios  Avanzados  en  Energla  y  Medio  Ambiente  (CEAEMA)  of 
Jaen  University.  Jaen  climate  is  characterized  by  a  high  level  of 
annual  direct  normal  irradiation  and  ambient  temperatures  that 
can  easily  reach  40  °C  in  summer  and  5  °C  in  winter  [59],  In  these 
climatic  conditions,  a  wide  range  of  cell  operating  temperatures 
was  registered  during  the  experiment,  what  allows  guaranteeing 
the  confidence  of  the  procedure  used  for  evaluating  the  accuracy 
of  the  different  methods.  Fig.  5  shows  the  histogram  of  measured 
cell  temperatures  during  the  two-years  campaign. 

From  the  collected  data,  the  parameters  of  the  different 
methods  were  determined.  Each  method  was  implemented  and 
the  root  mean  square  error  (RMSE)  in  C  between  calculated  and 
measured  cell  temperature  was  obtained  by 


RMSE(°C)=^(1/N)  £  (Tc 

,calculated  T c,measured  )2  (32) 

The  procedure  was  only  applied  to  methods  whose  outdoor 
validation  has  been  reported.  Results  are  shown  in  Table  3. 

As  can  be  seen,  methods  based  on  direct  measurements  on  the 
module  (RMSEs  from  1.7  to  2.5  °C)  are  more  accurate  than 
methods  only  based  on  atmospheric  parameters  (RMSEs  from 
3.2  to  4.3  °C)  as  expected.  Among  the  methods  based  on  direct 
measurements,  methods  based  on  module  heat-sink  temperature 
(RMSEs  from  1.7  to  2.2  °C)  have  the  same  order  of  error  magnitude 
than  methods  based  on  electrical  parameters  (RMSEs  from  1.8  to 
2.5  °C). 

Regarding  methods  based  on  module  heat-sink  temperature, 
the  method  of  1SFOC/IES-UPM  gave  lower  error  than  the  method 
of  Sandia  (1.7  °C  vs.  2.2  °C).  This  is  probably  because  a  detailed 
knowledge  of  the  materials  between  the  cell  and  the  heat-sink 
was  available  in  order  to  determine  the  “internal  thermal  resis¬ 
tance”  parameter  ( p )  of  the  1SFOC/IES-UPM  method,  while  the  AT 
parameter  of  the  Sandia  method  was  determined  from  an  average 
of  the  measurements  taken  around  the  reference  irradiance. 

Regarding  methods  based  on  electrical  parameters,  the  best 
behavior  was  observed  for  the  method  of  Peharz  (RMSE=1.8  °C). 
This  can  be  due  to  the  fact  that  the  functional  dependence  of  the 
cell  temperature  on  electrical  parameters  was  carefully  obtained 
by  the  authors  through  indoor  measurements  in  a  sun  simulator. 


Cell  temperature  range  (°C) 


Fig.  5.  Histogram  of  measured  cell  temperatures  during  the  experimental 
campaign. 


Table  2 

Characteristics  of  the  measured  HCPV  module. 


Geometric  concentration  Primary  optics 

Secondary  optics 

Type  of  solar  cells 

Number  of  solar  cells  Cooling 

500  x  PMMA  squared  flat  Fresnel  lens 

Refractive  truncated  pyramid 

Lattice-matched  GalnP/GalnAs/Ge 

6  Cells  in  series  Passive 
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Table  3 

Errors  between  predicted  and  measured  cell  tem¬ 
perature  for  different  methods. 

RMSE  (  C) 

Methods  based  on  module  heat-sink 
temperature 

Sandia  2.2 

ISFOC/1ES-UPM  1.7 

Methods  based  on  electrical  parameters 

IEC  60904-5  2.0 

Peharz  1 .8 

Fernandez  2.5 

Methods  based  on  atmospheric  parameters 

Almonacid  4.3 

Hornung  4.0 

Fernandez  ANN  3.2 


Method  of  IEC  60904-5,  which  is  based  on  physical  parameters, 
also  gave  good  results  (RMSE=2.0  C).  Method  of  Fernandez 
performed  slightly  worse  (RMSE=2.5  “C),  but  has  the  advantage 
that  it  is  based  on  a  simple  linear  expression. 

Regarding  methods  based  on  atmospheric  parameters,  the 
method  of  Almonacid,  which  is  based  on  linear  coefficients,  has 
the  highest  RMSE  (4.3  °C).  The  method  of  Hornung  performed 
slightly  better  ( RMSE = 4.0  °C)  probably  because  it  introduces  an 
exponential  correction  based  on  wind  speed.  The  best  of  these 
methods  is  the  ANN-based  method  (RMSE=3.2  °C),  which  is  able 
to  better  characterize  the  complex  relation  between  cell  tempera¬ 
ture  and  atmospheric  parameters. 


7.  Conclusions 

A  review  of  existing  methods  for  cell  temperature  calculation  in 
HCPV  modules  for  electrical  characterization  has  been  done. 
Methods  have  been  classified  into  methods  based  on  direct 
measurements  on  the  module  and  methods  based  on  atmospheric 
parameters.  Methods  based  on  direct  measurements  are  prefer¬ 
able  during  field  performance  analysis  while  methods  based  on 
atmospheric  parameters  are  interesting  when  direct  measure¬ 
ments  on  the  module  are  not  possible.  In  turn,  methods  based 
on  direct  measurements  have  been  classified  into  methods  based 
on  module  heat-sink  temperature  and  methods  based  on  electrical 
parameters. 

Five  of  the  reviewed  methods  present  the  advantage  that  their 
parameters  can  be  obtained  from  manufacturers'  data  or  from  the 
literature  (ISOFC/1ES-UPM,  Muller,  IEC  60904-5,  Ju  and  Yandt), 
what  makes  easier  their  application.  The  rest  of  the  methods 
require  an  outdoor  experimental  campaign  in  order  to  determine 
the  parameters  or  measurements  in  a  sun  simulator,  what  makes 
more  difficult  their  application.  We  have  also  identified  other 
technical  difficulties  related  to  the  application  of  some  methods, 
such  as  the  need  of  a  system  that  shutters  sunlight  (Muller),  the 
need  of  the  whole  I-V  curve  of  the  module  (Yandt)  or  the  need  of 
advanced  knowledge  in  ANN  (Fernandez  ANN). 

One  of  the  methods  is  intended  for  calculating  cell  temperature 
of  the  module  working  at  MPP  (Yandt).  Its  authors  predicted  a 
difference  of  up  to  14  °C  between  cell  temperature  at  open-circuit 
and  cell  temperature  at  MPP.  Thus,  investigating  the  cell  tempera¬ 
ture  at  MPP  could  be  useful  for  certain  studies  and  could  be 
interesting  in  future  research  works. 

A  comparative  analysis  of  the  methods  has  been  presented.  For 
this  purpose,  monitored  data  obtained  during  an  experimental 
campaign  carried  out  at  Jaen,  South  of  Spain,  on  an  HCPV  module 
has  been  used.  Results  show  that  the  accuracy  of  methods  based 


on  direct  measurements  (RMSEs  from  1.7  to  2.5  °C)  is  somewhat 
better  than  that  of  methods  based  on  atmospheric  parameters 
(RMSEs  from  3.2  to  4.3  °C).  However,  it  is  remarkable  that  the 
differences  are  not  very  large,  i.e.  methods  based  on  atmospheric 
parameters  are  a  useful  tool  when  direct  measurements  on  the 
module  are  not  possible. 

In  the  context  of  module  power  calculations,  the  obtained 
errors  for  all  methods  are  acceptable.  As  it  is  shown  in  Ref.  [40], 
if  we  take  into  account  the  less  accurate  method  (the  Almonacid's 
atmospheric  method  based  on  linear  coefficients,  with  a  RMSE  of 
4.3  °C),  the  error  induced  by  the  cell  temperature  estimation  in  the 
calculated  power  is  within  +  1.5%  of  the  real  power  in  88%  of  the 
measurements  carried  out  during  an  experimental  campaign  of 
two  years  at  Jaen.  This  means  that  the  less  accurate  method  can  be 
used  to  estimate  module  power  accurately  enough. 

However,  depending  on  the  specific  application,  some  methods 
can  be  more  suitable  than  others.  The  choice  of  the  most  suitable 
method  is  not  trivial  and  will  mainly  depend  on  the  availability  of 
module  information,  on  the  required  accuracy  and  on  technical 
issues. 
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